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ABSTRACT 



During recent years , logarithmically periodic antennas have been 
widely used due to their frequency response characteristics, simplicity 
of design and directivity. However, their theory of operation still is 
in a development phase, and very few models have been fully analyzed. 

The present paper is an attempt to analyze the operation of a zig-zag 
model that has the property of being symmetrical, and suitable for 
operation against ground. The radiation pattern of the antenna is 
obtained for different models of current distribution, and, finally, the 
impedance characteristics and an approximate current distribution are 
obtained, using non-uniform transmission line theory. The results 
obtained show reasonable agreement with experimental data, and confirm 
conclusions drawn from physical considerations. 
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I - LOG-PERIODIC ANTENNAS - HISTORICAL BACKGROUND 



1.1 - Frequency Independent Antennas and Log-Periodic Antennas 

Frequency independent antennas were originated in 1954, when V. H. 
Rumsey [14] started the idea that a structure characterized by angles, and 
without any characteristic length dimension, should behave independently 
of the frequency of operation. As a practical example, Rumsey proposed an 
equiangular spiral structure. J. D. Dyson carried out a study of such a 
structure composed of two plane sheets of metal extending from a central 
point according to this geometry, and he proved that the model not only 
had frequency independence characteristics, but also that the interesting 
part of the antenna was confined to the area inside an "active region," 
where the dimensions become of the order of half a wavelength. That is, 
the outward area is not necessary and the structure can be truncated at a 
convenient size, and will be useful within a range of frequencies whose 
lower limit is determined by the peripheral dimensions, and whose upper 
limit is determined by the spacing at the feeding terminals, near the 
center of the antenna. 

A structure like this is bidirectional and has circular polarization. 
Also, since the active region rotates along the spiral, the polarization 
rotates too. Strictly speaking, then, not all the characteristics of the 
antenna are frequency independent, but, instead, they repeat periodically 
every time the frequency is changed so as to reproduce a reference space 
condition. The input characteristics, pattern shape and field intensity 
are constant. 

One year later, in 1955, R. H. DuHamel, who was working with Rumsey 
in broadband antenna development, proposed a structure that, instead of 
being entirely defined by angles, was interrupted by discontinuities 
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Fig. 1 - DuHamel's log-periodic teeth structure 



8 



properly located- The first model built (Fig. 1) was formed by a 
circular sheet of metal, cut along two symmetrical circular sectors, and 
with teeth placed along the radial edges. The size of the teeth and 
their distances to the center, are related by a constant ratio p. It is 
evident that such a structure, within its size limits, will repeat its 
radiating characteristics every time the frequency changes by p; that is, 
for all the set [p^^f^] - A plot of these frequencies in a logarithmic 
scale will show them equally spaced by log p, and that is the origin of 
the name log-periodic given to this kind of geometry. 

The model shown in Fig. 1 is a bidirectional antenna with large 
half power beam width. Bending one of the halves of the antenna toward 
the other, so as to form an acute angle, D. E. Isbell obtained a unidirec- 
tional antenna with the surprising characteristic of radiating toward the 
vertex; that is, backward in relation to the feed. Meanwhile, proceeding 
with his work in the equiangular spiral antenna, Dyson devised another 
unidirectional antenna by wrapping the two metal sheets of the planar 
spiral around a cone. Again, he obtained a backward radiating structure, 
whose characteristics could be adjusted by means of the cone angle and 
the spiral rate of growth. 

These two pioneer models brought log-periodic antennas to the field 
of practical application; and they, and their variations, have been 
extensively used in the VHF band and above. They are particularly suited 
for military applications where more often than not, a single antenna is 
required to cover a very large range of frequencies. As it might be 
expected, these antennas have moderate gains, and for fixed circuit 
applications they do not match high gain structures. However, when a trade 
between even performance along a large band, and optimum performance only 
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at one frequency or narrow band is allowed, the log-periodic models offer 
the best possibilities. 

1.2 - Unidirectional Planar Structures 

Traditionally, it is with the lower frequencies that the search for 
efficient radiating systems has been less successful. Good electrical 
performance has been obtained at the expense of size, and, although 
theoretically it would be possible to build optimum systems, the costs 
involved in the ground occupied and the hardward employed make such attempts 
prohibitive. In the HF band, the best compromise between size and opera- 
tional characteristics, during the past three decades, has been the rhombic 
antenna. Today, the rhombic still is the practical structure with highest 
gain, and when used with a range of frequencies within a ratio of not more 
than 3:1, it still might be the best answer for an HF radio link. Its 
main inconveniences are the size and the variable take-off angle, but the 
latter can be used with advantage by appropriate selection of frequency 
according to ionospheric conditions. 

The log-periodic models of the type mentioned above are not, in 
general, feasible for HF long range communications. Besides being too 
bulky, they employ sheets of metal that would make the cost too high. 

Two important developments were to come, and solve these two problems. 

The first was when DuHamel and his research team verified that only 
the peripheral strips of the metal sheets are active in the radiation 
mechanism, so that, essentially, all the structures can be reduced to 
wire frames. 

The second was when Isbell [11] built the first array of log-periodic 
half-wave dipoles (Fig. 2a). Again, here a surprising fact was verified. 
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Fig, 2a - Isbell’s transverse dipole array 




Fig, 2b - A tapered-ladder for operation above ground 
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The structure, as all log-periodic structures so far realized, radiates 
backward; but, to obtain proper operation, every other dipole of the 
array is fed through a reversal of the feed line connections r 

This dipole array established the applicability of log-periodics in 
the HF band. The structure can be built for vertical or horizontal 
polarization, and in the latter case, by adequate choice of the slope of 
the feed line, either fixed or variable take-off angles can be obtained. 
The first models had tubular dipoles, with radius scaled to the structure 
period. Granger Associates, who might be properly recognized as the 
pioneers in large scale building of HF log-periodic antennas, have 
replaced the tubular dipoles by wire triangular elements with appreciable 
reduction in weight, making it possible to build the rotatable long range 
log-periodic. This is a steerable antenna, mounted on a pole high enough 
so that the take-off angles are as low as 10®, and light enough to be 
rotated by a small electric motor through a 360® azimuth span. 

Incidentally, one of the problems that arose with these antennas was 
the matching of the feeding line. The antenna terminals show a balanced 
input impedance of the order of 200 ohms, and, typically, a maximum VSWR 
with respect to 50 Q below 2:1 through the band of operation. Two hundred 
ohms is not the most convenient value to carry energy from a distant 
transmitter, and the problem was that of obtaining a matching device that 
would enable the use of 50 ohms standard coaxial cable. That need speeded 
up the research in the area of baluns and matching transformers with 
outstanding results. Granger Associates have supported their achievements 
in antenna production, supplying also baluns with the required broadband 
characteristics, and multicouplers that enable the simultaneous use of one 
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single antenna with a number of independent transmitters; the total power 
being the only limitation in the system plus a small separation in the 
frequencies of operation. 

Looking for more and better, the next step in log-periodics was to 
try to obtain a structure that would use a ground image with the purpose 
of reducing its size. The reversal in dipole connections excluded that 
possibility in Isbell's array, except by the use of inverting transformers 
between radiators, which does not seem to be a very practical idea. The 
more significant successful structures are the bent log-periodic zig-zag, 
the tapered ladder, and the monopole array. The bent zig-zag, developed 
by J. Greiser and P. Mayes [8], [9], will be the main subject of this 
paper. The tapered ladder, developed by Wickersham, consists of an array 
of grounded vertical radiators (Fig. 2b) capacitively coupled to an 
inclined transmission line, the slope of which is adjusted for proper feed 
of each radiator. A variation of this arrangement, that uses triangular 
wire radiators, has been commercialized by Granger Associates. The 
monopole array, developed by D. Berry and F. R. Ore [23] uses phasing 
stubs between every two radiators, scaled with the structure period. The 
stubs are horizontal and open-circuited, so that only one wire is actually 
necessary, since its image provides the other, and the radiators are 
vertical tubes, of a quarter of a wavelength, for the respective operating 
frequency. 

In general, the half-size structures do not have as good a performance 
as the transverse dipole, and the recent advances in construction 
techniques have made the latter a stronger competitor, even for frequencies 
as low as 2 MHz. 
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1.3 - Developments in the Theory of Log-Periodic Antennas 



The development of log-periodic structures during their first years 
was mainly the result of good intuition and extensive laboratory testing 
of a large number of proposed models • As a rule, only after a successful 
model was obtained, was a theoretical foundation to explain its perform- 
ance developed, and only relatively recently models based on sound 
theoretical investigation are being realized. 

There have been three main approaches to the problem of analysis of 
operation of LP’s. The first, which we can call the classical approach, 
uses the theories developed for the components of the antennas and combines 
them to suit the particular model under consideration. The outstanding 
example of this approach is Carrell’s analysis of the transverse dipole 
array [10] , who described the behavior of the structure in terms of 
currents, voltages, input impedance and far-field with very close agree- 
ment with experiments. This approach usually has the limitation of being 
valid for a particular structure only. 

Another approach, that has been pursued by DuHamel [15], [16], 

Armstrong [16], Mittra, Jones [17], [18], [19], and others is to develop 
an adequate body of information and knowledge in log-periodic networks, 
on which a general theory of the antenna problem could be based. Great 
progress has been made toward this purpose, and results obtained have 
helped in a better understanding of log-periodics , and, in some cases, 
provided new ideas for future realizations. 

The third approach, followed by Mayes, Deschamps, Patton [20], Greiser 
[8], [9], Oliner, Mittra, Jones [22] and others is an application of 
concepts developed for surface and leaky wave structures, that are 
typically uniformly periodic structures, to the active region of log-periodics 
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that can be approximated as a uniformly periodic region for limited 
purposes. Since such concepts give an excellent insight into the main 
operational characteristics of these structures, we will briefly review 
some of its basic points [4], [5], [7]. 

Consider a structure, such as a waveguide, with a periodic grating 
in one face, and suppose that the waveguide carries an electromagnetic 
wave. Due to the discontinuities in the grated face there will be 
leakage of energy through them, and, also, the fields inside the guide 
will adjust to this periodic configuration. According to Floquet’s 
theorem, ”for a given mode of propagation at a given steady-state 
frequency, the fields at one cross section differ from those one period 
away only by a complex constant" [7]. This can be expressed as 



E 



2 






-YL 



where L is the period of the grating and y is the propagation factor. 

Now, if the longitudinal dimension of the guide is made to coincide 
with the z axis, and we expand the equation in a Fourier series, we get 



E = E(x,y,z)e"^^ = I E^(x,y)e"^ 



( 1 ) 



all n 



that is, we obtain a set of n space harmonics, each one having propagation 
constant 



' 27Tn 

r = Y + j — . 



If we now assume that the system is lossless, we obtain 






( 2 ) 



where 6^ is the medium phase constant. Now, we apply this result to the 
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usual description of waves bounded by guides. In general, the fields 
are described by the wave numbers and a relation of the form 



2 2 2 2 
K = K + K + K . 
X y z 



When the x = 0 plane is a ground plane, that is the general case 
of interest for the present purpose, the y number is zero, and, with the 
assumption of no losses, we have 



6 



2 



K 




+ 6 



2 

z 



where 6 = K is the free space phase constant. 

If now we recognize that the x and z phase constants are given, for 
the grated guide, by a series of the form of equation (2), we obtain the 
general form of the solution. But, before, let us consider an angle 0, 
as being the angle between the z axis and the direction of a z component 
wavefront. Then we can relate the x and z wave numbers, or the phase 
factors as 



6 =» K sin e 

X 

and 

6^ = K cos 0 , 

so that the direction of leaking energy, that is the direction of the 
wavefront, is given by 

6 

cos 9 = • (3) 
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Now, according to (2) we recall that 3 can take any discrete value 

z 

6^ (note that n ranges from all negative values to all positive values), 
and, being so, the ratio (3) can take an equal number of values, positive 
or negative. When this ratio falls in the range -1 to +1, that corresponds 
to a solution of the cos function, and also to real solutions of space 
liariuonics for which there is leaking energy. In particular, when (3) takes 
the value -1, energy is fired backward, and this corresponds to the mode 
of interest in the present analysis. 

Since the amplitude of the harmonics decreases with increasing order, 
we would like to have the structure operating in a low order mode. For 
the n = -1 mode we have 



and 



cos 




^0 - 2tt/L 
K 



In log-periodic structures, always there is some sort of guide that 
carries energy from the feed point toward the radiators, and discontinuities 
placed along the guide force energy to leak away. In the transverse dipole 
the line is the feeder, a transmission line, and the grating is represented 
by the dipoles,* in the teeth structure of Fig. 1, the plane sheets form 
the guide and the teeth form the grating; the bent zig-zag is a transmission 
line deformed by the radiating elements. In each of these models, and, in 
general, in each log-periodic model, we can consider three distinct 
regions: the active region comprising the two to five radiating elements 

near resonance for the operating frequency, the feeding region, between 
this and the source where radiation is small and in general the radiators 
can be considered as a capacitive loading, and the end region, that we 
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assume is not reached by any significant amount of energy, since most 
of it was radiated by the active region. Since the active region is 
small, we now assume that it is a uniformly periodic structure, and then 
apply it to the above results for leaky waves. 

Perhaps a clearer picture of such an application is by means of 
dispersion diagrams. Fig. 3 shows a dispersion diagram for equation (3), 
assuming that Sq/K >1. It is a well-known fact in transmission theory 
that propagation within boundaries takes place when the phase constant 
in the medium is larger than the phase constant in the boundary material, 
i.e., when the propagation velocity in the medium is lower than the 
propagation velocity in the boundary material. For metallic boundaries 
the velocity is usually taken as that of free space, and we represent 
that condition by the phase constant K. This situation is usually 
identified by the designation of slow waves and it is said that periodic 
structures have stop bands for fast waves. These fast waves that do not 
propagate within the structure are the ones that leak away and form the 
radiated field. 

In Fig. 3 we represent the fundamental and the first and second 
negative harmonics as given by equation (2) , and we show the +45® slope 
lines that separate zones of slow waves from zones of fast waves. We 
see that the fundamental is a slow wave for all frequencies, and so are 
all the positive order harmonics. The negative order harmonics are 
slow waves for low frequencies, but as the frequency increases to near 
the point B, the first negative harmonic approaches the negative slope 
line. On the line we have 

8q-2it/L 
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Fig. 3 - Dispersion diagram 
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that is, we start having radiation backwards; and until the positive 
slope line is reached the radiation angle rotates from 180° to 0°, 
Meanwhile, the other negative order harmonics were having similar 
behavior in the region that we can visualize to the left of the diagram. 

Since the analysis is applied to a uniformly periodic structure, 
we only have well-defined active bands, the first of which is the 
strongest, separated by stop bands (these have the opposite meaning of 
the stop bands of transmission theory). However, in log-periodic 
antennas the structure is scaled almost continuously (typical scale 
factors range from .80 to .98), so that there is a continuous active 
band, limited only by the operating frequencies of the first and the last 
of its elements. The frequency independence characteristics of the 
antenna depend on how well its discrete scaling approaches a continuous 
scaling, and the over-all antenna performance can be predicted from this 
type of analysis. Greiser and Mayes [8], [9] have used this technique 
with remarkable success, in predicting whether or not different models of 
bent zig-zag antennas would radiate, and how the radiation characteristics 
could be improved. Mittra [17] has conducted a similar analysis for the 
transverse dipole, showing by this method how the phase reversal was 
necessary to make the dispersion diagrams approach the radiation zone. 

Again, there is a remarkable lack of dispersion diagrams developed 
for the type of structures used in log-periodics , and, as yet, the theory 
was not really adapted for them, except under the heuristic approach 
described above. We should note here that, although this application of 
leaky waves theory in a wired log-periodic structure might, at first 
sight, seem rather displaced, there is reason to believe that the theory 
is actually valid. R. L. Bell, C. T. Elfving and R. E. Franks [21] 
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reported measurements of fields around a log-periodic structure that 
show the existence of a forward wave progressing from the source along 
the structure, and a backward wave with a 90 degrees special phase 
difference. The forward wave is supported by the structure; the backward 
wave is launched by it into space. That is exactly the type of leaky 
wave behavior used in the above theory. 

The mathematical tools developed to deal with the various problems 
posed by electromagnetic theory, although having in common the basic 
source given by Maxwell’s equations, follow different courses and tend 
to place in distinct categories, problems that' have much in common. It 
seems to be a marked trend in recent years, to reunite electromagnetic 
theory in the sense of extending concepts and ideas developed for 
particular systems to other, apparently different, applications. 

As a final remark, we note that, since log-periodics are scaled up 
from the feeder, the active elements are the ones in the active region 
together with the smaller ones. Smaller elements in all cases tend to 
increase the propagation velocity (the feeder approaches a uniform line 
toward the feed point). If we would attempt to build a log-periodic for 
forward radiation we would scale the structure in the opposite sense, 
and would count on positive order waves for the radiation process. 
However, while in the usual scaling the positive harmonics find a high 
impedance in the first part of the feeder (they are in the stop band of 
transmission theory and mismatched to the available radiators), in the 
opposite scaling the negative harmonics would still reach the end part of 
the feeder and cause a disturbing radiation from that region. 
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II - THE BENT LOG-PERIODIC ZIG-ZAG ANTENNA (BLPZZ) 



The BLPZZ (Fig. A-1) was developed by Greiser and Mayes using the 
concepts of dispersion diagrams. As mentioned before, this was one of 
the first attempts to obtain a LP antenna with S 3 mimetrical characteristics 
that could be fed against ground. For that reason, it seems particularly 
suited for HF applications. 

The basic geometry of the antenna is derived from a zig-zag antenna 
bent along a convenient line, so as to obtain adequate radiating elements 
and proper connecting stubs. Since the structure is progressively 
elevated from ground the connecting stubs are actually radiating, but, 
for small elevation angles that can be neglected. Their essential function 
is to provide the proper phase delay between radiators, and they can be, 
and have been, replaced by other forms of delay line, or even removed, 
without destroying the basic characteristics of operation of the antenna. 

The antenna is usually defined by the radiators spreading angle a_, 

£i 

the stubs spreading angle a , the elevation angle e, the period p, and 

s 

the length of a wire, 2b, or the distance of one element to the geometric 
origin. In Appendix A we develop the relations among these parameters and 
other useful quantities. We note that the spreading angles and the period 
are interrelated through the slope of the wires. Although each one of 
them can be independently changed, that involves a change in the dimensions 
and inclination of the wires. Also, we note that each element, for small 
elevation angles, can be compared to a rhombic antenna made of small 
wires, and with a very large apex angle. 

This antenna has the interesting characteristic of being fed in 
series, i.e., the radiators are part of the feed line. The other planar 
LP structure so far analyzed — the transverse dipole — has the radiators 
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fed in parallel, and it might be worthwhile to compare the results 
obtained. On the other side, this one, also is a structure that seems 
adequate to analyze with conventional tools. 

In the first part of this work we obtain a set of equations for the 
radiation pattern of the antenna. Then we assign to the wires different 
types of current distribution, and compare the results with the 
experimental data obtained by Greiser and Mayes. 

In the last part of the work we consider the antenna and its image 
as a non-uniform transmission line and we obtain an approximation for 
the input impedance and current distribution. 

2*1 - The Radiation Pattern Equations of the BLPZZ 

To obtain the radiation pattern of the BLPZZ structure above ground 
we begin by considering one element at a time. We define as an element, 
one period of the structure. Since the elements are related by a constant 
size factor, that is, the period p, the general equations obtained for 
one element are valid for the other elements as well. 

In each element, we assume that the connecting stub does not radiate, 
and we split the radiating section in two straight wires. Then we 
obtain the fields due to these two wires and their images, all referred 
to an element phase center, and add them vector ially. 

We make two more assumptions. First, we assume that there is no 
mutual coupling between elements; second, we assume that the wires carry 
travelling waves. This corresponds to considering the structure and its 
image as forming a radiating transmission line properly terminated; that 
is, without standing waves. Later we will further elaborate on this 
point. 
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After obtaining the element’s equations we establish an array 



expression, accounting for spacings and feed currents. 

We outline here the development of the equations for one wire, wire 
1 in Fig. 4. Details, as well as the work concerning other wires, can 
be found in Appendix A. 

In general, the distant field due to the flow of current along a 
wire is given by 



where 0’, s, 1 and r are defined in Fig. 4, and w, k, and it have the 
usual significance. p(s) is the current distribution. In general, this 
will be given by 



where a is the attenuation factor and 3 is the propagation factor, along 
the wire. 

For large distances the angle 0 and the distance r are essentially 
independent of the position along the wire, so that we can consider them 
as constants in regard to s, provided that we account for the phase delay 
at the far field due to the different points along the wire. 

Again, referring to Fig. 4, if we consider the distance R from the 
center of the wire, we can express r as 




(4) 



m(s) = Iq e 



-(a+j3)s 



(5) 



r = R - s cos 0 



and, replacing in equation (4), and defining 




we obtain 
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Fig. 5 - Geometrical relations in the BLPZZ 
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Now we apply this equation to wire 1 in Fig« 5., Here 1 = 2b, and, then 
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Defining 



wPq ^-j(kR+wt) 
Yk R 



-j (kR+wt) 

60 I 



we obtain 




[«nKl 

e - e 

[®< t- ^ k(~f‘ ■“ ^i)l 




Since we are going to use this equation with a digital computer, this 
is a convenient form*, 

Now we need to convert this equation to a system of coordinates to 
be used with all the structure* We choose spherical coordinates, as 
shown in Fig* 6, and we obtain the following transformation equations (see 
Appendix A, Parts II and III): 
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( 8 ) 



COS 0^ = cos (ye) cos 0 - sin (y-c) sin 0 cos 0 



E 



01 




(9) 



_ p sin (y-c) sin 0 
'0^ ■ ^0| sin 0| 



( 10 ) 



At this point it is convenient to introduce the correction to the 
"phase centers" of the element. Here the designation "phase center" is 
used to mean the reference point of the element. The actual phase center 
was not computed, although physical considerations would lead to the 
point C in Fig. 5, that was chosen. 

Again, we assume that the angles 0 and 0, as well as the distance R 
are not altered to an appreciable degree, and we only take into account 
the phase delay in the far field. The distance c^, from the center of 
the wire to point C is given by (see Appendix A, Part I) 



^ b sin (y-£) 
1 cos 6^ 



( 11 ) 



where 



6 



1 



-1 L sin £ -f b cos (y-s) 
b sin (y-e) 



The phase difference in the far field is given by (see Fig. 5) 



- R = c^ cos a, 

where 



cos 



[sin 6^ cos 0 + cos 6^ sin 0 cos 0] , 
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Fig. 7 - Correction for reference point 
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so that the factor e in equation (6) becomes becomes e ^ 

and, finally, combining equations (7), (8), (9), (10), and this correc- 
tion, we obtain 




-f-C,Co5o<, 



) 



C — 0 



X 











and 




-hC|C0Sc<j 



) 



£. - e 

( 15 ) 



Now combining these equations with the similar equations obtained 
for the other three wires (see Appendix B) we have, for each element. 



"e = "e, + "e, "e, + \ 

n 1 2 3 4 



(14a) 



»„ "l h ®3 ®4 



(14b) 



For the whole structure we will have 



N 

"e- I Ve 

n=l n 



(15a) 



N 

= y a E. 

0 n 0 

n=l n 



(15b) 
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where is a factor accounting for the spacing and feed currents, that 
can be split in the product of two terms, one accounting for each. 

The feed current for the n+1 element is the feed current for the 
order n element delay attenuated by the length of wires 1 and 2 and the 
connecting stub. As mentioned before, the connecting stub is assumed 
to be a uniform transmission line, so that it is only responsible for an 
additional delay. 

Referring to Fig. 5, the distance along the wires from the feed 
point of element 1 to the feed point of element 2 is 






From Appendix A we have 

b 



_ _1 tan as 
1 /-* tan aE 

/p 



so that the distance becomes 



dj . 4bi 1 + 



tan as 
/p tan aE 



The distance to the feed point of element 3 is 



S = 



1 + 



tan as 
tan aE 



+ 4b. 



1 + 



tan as 
/p tan aE 



4b, 



1 + 



tan as 
tan aE 



1 + — j , since b^ = b^^; 



and the distance to the order n element is 



d = 4b, 
n 1 



1 + 



tan as 
/p tan aE 



1 ^ i . ^ + 1 

p 2 n-2 

p p 



= 4b, 



1 + tan as 



SUM 



/p tan aE 
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where SUM = 1 + — + + . . . H r 

p 2 n-2 

P P 



If we subtrace — from S, where, 
P 



p p 2 3 n-l’ 

P P P 



we get 




_ 1 . s = - 

n-1 ^ ^ 1 - 1/p 

P 



and the distance becomes 



d = 4b 
n 



tan as 
/p tan aE 



1 - 1/."-^ 

1 - i/p 



If we call 3 the propagation constant along the radiating wires, and 
assume that along the stubs the constant is K, as in free space, we have 
the phase angle of the feeding current given by 



4b^K 




tan as 
/p tan aE 



1 - i/p^'^. 

1 - i/p 



(16) 



We could obtain an analogous expression for the attenuation, assuming 
that the attenuation factor was constant along the radiating structure; 
however, physical considerations imply different attenuations along the 
line, so that we obtain, for the order n element feed current an attenua- 
tion given by 

n-1 

y 4b. a. (17) 

i=i ^ ^ 



where is the average attenuation for each element. 

Then the correction factor for current in the array expression is 
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Fig. 8 - Correction for array spacing 
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e 




(IS) 



To introduce now a correction for spacing, we are going to use as 
the reference point the '‘phase center" of the array. Again, we do not 
know exactly where the phase center is, but we expect in that respect a 
behavior similar to other log periodic structures [10], [21]. Being so, 
the phase center should lie by the element that is nearest to resonance, 
and that we call the M element. Then, according to Fig. 8, we obtain 



= R - (L^ cos e - cose) cos 0" 



where cos 0" = sin 0 cos 0 and 



R = R - “ L ) sin 0 cos 0 cos e 

n M n 



(19) 



^ -ikR , +ik(L. “L )sin 0 cos 0 cos e 

In equation (6) , e becomes e e ^ n 



and we insert in the array expression the factor 



e 




( 20 ) 



Then the total correction a is 

n 




C^O 



Finally, the array equations become 
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>>=1 



(2 2a) 



^0~^, ^ 






cos £. 



’< £;a» 

(2 2 b) 



And from here we obtain the total field E by taking the real parts of E 

0 

and E^ and combining them vectorially, so that 



E (&,0) = 




Reati^ejf +[Reat(^0)f 




Real (^d ) 
Red(fei) 



= \BlH 



where ¥ is an angle counted positive in the clockwise sense starting in 
the direction of the vector by an observer at the array phase center, 

2.2 - Current Distribution - Radiation Diagrams 

Before programming the equations obtained in section 2.1 in a 
computer, we had to decide on the current distribution that would be 
most closely approximate to the real situation. With that purpose in 
mind, we made the following considerations: 
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Note: n denotes resonant element 



Fig. 9 - Models of current distribution 

(1) As obtained in other LP networks 

(2) Assuming travelling wave attenuated at 
active region 
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(1) The antenna and its image form a radiating transmission 
line, and since experiments [8], [9] show that its input impedance 
is mainly resistive and almost constant along the operating band, 
having a small capacitive reactance probably due to input termination 
effect, it is reasonable to assume that the line is well terminated, 

(2) The elements forming the feed region must see the same 
impedance toward the end, and that is the input impedance of the 
antenna. This is so because that impedance does not essentially 
change with frequency, and a change in frequency is equivalent to 

a change of antenna. That is, when we decrease the frequency by one 
period, we move the active region one period toward the end, and 
that is equivalent to using one antenna ’*one period larger;" the 
impedance that was seen by the input terminals must now be seen at 
the end of the first element. But since the input impedance remains 
the same, that means the first element is a good match, and so must 
be all the elements forming the feed region. Then the current along 
the feed region must be essentially constant, except for a relatively 
small amount of attenuation. 

(3) The energy reaching the end region must be very small. The 
structure is abruptly terminated without any dissipation load, and 
it has been shown experimentally [8], [9] that the effect on the 
radiation pattern of ending the structure by a short or an open 
circuit is very small, for frequencies well above the lowest frequency 
of operation. Then, it seems reasonable to assume that the current 
dies off after the active zone, and is essentially zero at the last 
elements . 
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(4) The current distribution at the active zone is more 
problematic. The non-uniform line forming the antenna is continuous 
in the sense that its impedance [2] changes without interruption, 
and has the property of a constant product of reactance times 
susceptance. We can call this a continuously tapered line, with 
linear taper, of the type often used to match impedances. Under 
this point of view it is continuously matching itself. 

On the other side, if we think of the line as composed of series 
sections connected by stubs, then we have to consider the fact that 
some of these sections (the ones forming the active region) are near 
their natural resonance frequency and are likely to form standing waves. 
The first point of view would lead us to a heavily attenuated 
travelling wave. The second leads us to the type of current distribu- 
tion found by Carrel [10] in the dipole array, and by Mittra and 
Jones [18] for a theoretical line, A point to note here is that both 
these cases had parallel loads, while the present case is of series 
loading. Fig. 9 shows the two types of current distribution corre- 
sponding to the two points of view. 

We started by trying these two models. The results obtained can be 
summarized with reference to Fig. 10, 11, 12, that show diagrams obtained 
with the data given for the antenna referred to by BLPZZ 19c in Greiser’s 
and Maye’s work. This antenna operates from 500 to 2500 MHZ. Fig. 10, 
for 2000 Mcs is what can be called the typical pattern. Most computer 
runs show that kind of diagram. Half power beam widths of the order of 90® 
in azimuth and 35® in elevation. The directivity, computed by the 
approximate formula given by Kraus [6], is around 7,5 db. Also, in almost 
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every case the difference between the results for the two types of 
current distribution is small, although the distribution of curve 1 in 
Figc 9 gives slightly better results, particularly in what concerns the 
f ront-to-back ration 

After an extensive set of runs was obtained, it was found that for 
certain frequencies the diagrams depart appreciably from what was 
expected. Two extreme examples are shown in Fig. 11 and 12. The first 
is ”too good,” and the second is ”too bad,” The most likely explanation 
that we can find for that is in the way the current distribution is 
processed in the program^ The program has a sample of the curves in 
Fig. 9 and assigns a value to each element. Now, that sample is always 
the same, but the real situation in the antenna varies with the frequency. 
When it happens that the wires forming one element have length of one 
quarter of a wavelength the sample fits exactly, but when that dimension 
”falls” somewhere between two elements, it is clear that we are doing a 
rough approximation, because in that case the program takes the element 
just below it, and treats that element as if it were a quarter of a wave- 
length. 

To explore this point further, we made a set of runs at closely 
spaced frequencies, covering one period of frequency values. Table I 
shows the results. For each frequency, column 2 shows the corresponding 
quarter of a wavelength, columns 3 and 4 show the sizes of the wires 
nearest to that value, columns 5 and 7 show the obtained directivity and 
orientation of the E vector all around the ground plane for the current 
distribution of curve 1 in Fig. 9, and columns 6 and 8 show the same 
results for the current given by curve 2. 
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3 


4 


5 


1 

6 1 

i 


7 


8 


Freq 

MHz 


X/4 


Wire’ s 


length 


Current Dist 
(1) (2) 1 
DIR/ISO DIR/ISO! 


Polarization 
(1) (2) 


820 


.091 


.078 


.092 


5.24 


5.00 


-90° 


-90° 


855 


.088 


.078 


.092 


1.45 


7.05 


+90° 


+90° 


875 


.086 


.078 


.092 


6.36 


6.06 


+90° 


+90° 


895 


.084 


.078 


.092 


5.24 


3.9 


+90° 


+90° 


935 


.080 


.078 


.092 


8.89 


10.9 


+90°within 

+60° 


-90° within 
+60° 














+90°remain. 

directions 


+90° remain, 
directions 


970 


.077 


.067 


.078 


7.05 


7.45 


-90°within 

+135° 


-90°within 

+130° 














+90° remain, 
directions 


+90° remain, 
directions 



Table I 

Variations in Results Along One Period of Frequency Values 
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Here it becomes clearer that the current distribution 1 results in 



a more balanced set of results than distribution 2. It is interesting 
to see the way the direction of the field vector changes along the period. 
Actually, this change of frequency is equivalent to moving the structure 
one period in space, at fixed frequency, and, in any case, we would 
expect the change in phase that has been observed in all LP’s. It has 
been called the phase rotation principle — the phase of the resultant field 
changes by 360° with a change of one period, in space or frequency. 
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2.3 The BLPZZ Antenna as a Non-Uniform Transmission Line 



With the purpose of obtaining a better description of the current 
distribution in the BLPZZ, we consider the structure and its image as 
forming a non-uniform transmission line. For easier handling of the 
equations involved, we split the structure in sections, each composed of 
a radiator element and the following connecting stub. Then we apply to 
the radiator Schelkunof f ’ s equations for slightly non-uniform lines, 
and assume that the stubs are uniform lines with characteristic impedance 
equal to their average impedance. To account for radiation, we introduce 
in the equations an attenuation factor and a phase factor obtained from 
the radiation resistance of each element and its average impedance. The 
value of radiation resistance used is one obtained by Chaney [24] and Klamm 
[25], for rhombic antennas in free space. 

According to Schelkunoff [2] a non-uniform transmission line can be 
described by equations of the form 

V(x) = Vq(x) + V^(x) + V^Cx) + ... (25a) 

I(x) = Iq(x) + I^(x) + I^Cx) + ... (25b) 

where Vq(x) and Iq(x) are the expression for a uniform line with the 
average characteristics of the line under consideration. When a line or 
a section of line is slightly non-uniform the terms of order 2 and above 
can be neglected, without large error. In particular, when the line has 
the property of having a constant ZY product, where Z is the line 
reactance and Y is the line susceptance, and Z and Y vary slowly about 
their average values, the first order term is given by 
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( 26 a) 



Vfl*) iVIjLBw coshT X s'lTih Px] - Kp Ij,[/\(x)cos^rx-B(x) smAPx] 



Ii (x')ri^ [jB(x)sm<)rx-A(^)cosKPx] t siriKr)<-B(^)tos^rx] (26b) 

Ko 



where 



r = /zY 



is the propagation factor 



is the average impedance of the line 
A(x) = — / Z (s) cosh 2Fs ds 

^0 0 

, X ^ 

B(x) = — J Z (s) sinh ZTs ds 
0 0 

Z = Z - Z 

av 

Z = Z is the average reactance of the line 
u av 

Yq = is the average susceptance of the line, 

and we assume that the wires are perfect conductors, perfectly isolated, 
so that there is no resistance or conductance to account for. 

The required condition of slow variation of line parameters can be 
met by considering conveniently small sections of line at a time. 

Using the hyperbolic form of the equations for a uniform line, we have 

Vq(x) = Vq cosh Tx - sinh Tx (27a) 

I (x) = - 77 “ sinh Fx + cosh Fx (27b) 

0 U 
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and, then, 



= V, [( I B (x)) cosh r< - A r'x ] - Ko lo CO S ll fx + (l - B(X )) si TO K Px j 

{ZB a) 



J(X) = ::^[(l-t-B(x))5mArx-A(x)cosAr'xJ+Jo [A(x) sinliPK 4- (l-B(x)) Cosh Px j 

^ (,2 2 k ) 

For the line parameters we use the equations for the inductance 
and capacitance per unit length of a two wire line, in the logarithmic 
form, implying the assumption that the radius of the wires is very small 
compared with their spacing, and, so, 

Z = j^ln(-) = i 120k In ( - ) (29) 

7T a a 

Y = j w e 7r/ln ( *^ ) = j k/120 In ( ~ ) (30) 

a a 

where 

w = 27 t times the frequency 

y = 47T X 10 ^ Henry/m is permeability of free space 
-9 

e = 10 /367T Farad/m is permittivity of free space 

s is the spacing between wires, and 
a is the wire radius. 

The propagation factor, F, is a constant if we do not account for 
losses. If we consider the radiation resistance as a distributed loss 
resistance, with value Ra per unit length, we can use the following 
approximations described in the literature for lossy transmission lines, 
and define an attenuation per unit length. 



a 




( 31 ) 



and an average phase factor 
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6 = 



(32) 



/ 1 1 


1 1 


M- + ~T 


= k 1 + -^ 


' 8Q^ 1 


' 8Q^ / 



where 

R 

^0 

k is the free space phase constant, and 
r = a + j6 

As mentioned before, we use for the computation of radiation resist- 
ance the equations given by Chaney and Klamm for a rhombic antenna in 
free space. Those equations were developed for an antenna well terminated, 
and carrying a travelling wave. Actually, the present structure does not 
exactly have the rhombic shape, but the elevation angles considered are 
small enough to allow the assumption that they can be neglected. More 
daring is the assumption of a travelling wave. However, from the work 
reported by DuHamel on nonuniform transmission lines [16], we conclude 
that the value of Q should not be expected to be critical, and that is 
essentially what would change directly with the value of radiation resist- 
ance. 

So far we have not accounted for the mutual impedance among radiators, 
and, naturally, that is due to the large amount of computations needed, 
that would exceed the time allowed for the present work. After presenting 
the results obtained with the present assumptions we will try to do a 
rough estimate of the mutuals effects, based on the values found for 
current distribution and impedances. 

Using equations (28) for the radiators and equations (27) for the 
stubs, we start from the end of the line, where there is a short-circuit, 
and proceed toward the feed point computing the input impedances, and the 
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coefficients of the equations. As a block the equations can be expressed 
by 




(33a) 






(33b) 



and, for each section, we obtain A, B, C and D, and the ratio 

When we reach the feed end of the line ratio V^/I^ is its input impedance. 



Then we assign to the current the value 1.0 + jO.O and we proceed forward, 
using the previously computed values of A, B, C and D, obtaining the 
values of current and voltage along the line, for as many points as the 
number of sections considered. From these, and using the equations, the 
complete description of the line can be obtained. 

We program these equations into the IBM 360 system, according to the 
program in Appendix C. Appendix D shows the details of the equations 
involved. 

Sample program SMIGUEL shows the results for the same antenna that 
was used in sample program AZORES (radiation pattern) , for three frequencies 
within the operation band — 1100, 1700 and 2300 MHZ. For each frequency 
the first run does not account for mutual impedance effects. The second 
run includes the mutual resistance that would be found between parallel 
wires with sinusoidal current distribution. 

The experimental results obtained by Greiser and Mayes [8], [9], with 
the BLPZZ fed against ground show that the input impedance has a resistive 
component varying between 90 and 250 ohms and a reactive component 
varying between +10 and -100 ohms. Since we are considering a double 
structure, we would expect the resistance values to be about doubled^ 
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Although the behavior of the reactive component is more difficult to 
predict, we would expect the same orders of magnitude. 

These assertions are confirmed in the results obtained for 2300 MHz. 
Referring to SMIGUEL the column of ZINR values shows the input impedance 
seen at the feed end of each radiator, and, not only the input impedance 
is as expected (291. 3-j 314 . 8 ohms), but also, as assumed before, each 
element essentially provides a match for the next element. The column of 
CUR values shows the magnitude of the current at the feed end, at the 
vertex, and at the stub end of each radiator. We note that, for this 
frequency, element 2 is smaller than a quarter of a wavelength, and that 
element 3 is larger than that. Up to the third element the current 
magnitude values are typical of an attenuated standing wave situation. 
After element 3 the distribution is of a heavily attenuated travelling 
wave. 

As mentioned beofre, the more delicate assumption in this approach 
was to apply to this structure the radiation resistance given by Chaney 
and Klamm [24], [25] for the rhombic antenna. It appears from these 
first results that such assumption is reasonable for the large elements, 
but it is not adequate for the small elements. Also, for these elements, 
the effects of mutual impedance become increasingly important, since 
they are "electrically closer" to each other. These two sources of error 
become more important as we lower the frequency. At 1700 MHz (the 
resonance dimension falling between elements 7 and 8) the results 
obtained for the first 4 elements are inconsistent. 

A closer inspection of the results shows that the departures in 
impedance values are related to extremely high values of Q (above 20) , and 
that fact even more strengthens the idea that the loading due to mutual 
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effects plays an important part in the behavior of the "feed region" of 
the antenna. 

At this point we found it useful to make some numerical analysis of 
radiation and mutual resistance values, with the following conclusions: 

(1) The radiation resistance of the "rhombic" elements increases 
steadily with size, from a few ohms at the "feed region" to 
thousands of ohms at the "end region." 

(2) The mutual resistance between parallel wires in adjacent 
elements with sinusoidal distribution of current, in phase, decreases 
steadily from hundreds of ohms at the "feed region" to essentially 
zero at the "end region." 

This enables us to obtain a better approximation, by also including 
the mutual effects, taking two parallel consecutive wires of the structure 
at a time. Again here, we neglect a few points; namely, that the parallel 
wires considered do not have the same size, the input currents are not in 
phase, the current distribution is not exactly sinusoidal, and that there 
are other interactions to be accounted for, among the non-parallel wires. 
However, before taking this step we made extensive numerical calculations 
with different arrangements of wires with different conditions of phase 
differences and current distributions, and the results obtained were all 
within the same orders of magnitude. 

Essentially, what we are doing here is to introduce a rough correction 
to lower the Q of the "feed region," and the same result could be obtained 
if we assigned to it specific (and reasonable) values. As shown in the 
referenced paper by DuHamel [16] the exact value of Q is relatively 
unimportant. 
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We note that, by this procedure 

(1) at the "feed region" we are essentially neglecting the 
travelling wave approach since the rhombic radiation resistance is 
much smaller than the mutual resistance, and 

(2) inversely, at the "end region," we are neglecting the sinus- 
oidal distribution approach* 

The results obtained with this procedure apparently prove its 
soundness* For the three cases previously considered, the input 
impedance values now fall in the expected range (439. l-jl48* 5, 562.1-jl09.8, 
and 516*4+jl9.9 ohms)* The current distribution shows the same features as 
before but with stronger attentuation, making it resemble more closely an 
attenuated travelling wave configuration. 

The next step would be to obtain a more exact description of the line 
by considering a larger number of smaller sections, obtaining better 
approximations for the non-uniformity corrections, and a more detailed 
picture of the current values. That would enable the application of these 
results to the radiation equations. Unfortunately the attempts to program 
this detailed analysis failed due to unknown reasons, related to the 
difficulties encountered during the "adaptation" phase of the new IBM 360 
system, recently installed at the Naval Postgraduate School. Although we 
have received extensive help from the Computer Facility staff, and 
undoubtedly the program would eventaully work, there was not time available 
to pursue the project* 

Before closing this subject we note that, due to the method used in 
deriving the coefficients in equations (33), there is a cumulative source 
of error in the approximations used* This is so because every new value 
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is based on another value that had been obtained by an approximation, 
and only the first value is an exact value. Since we start from the end 
of the structure, the ”end region** coefficients are little affected, but 
the ones at the **feed region** are likely to be farther away from exactness. 
This might be another reason for the worse results obtained with the 
smaller elements. The way to correct, or at least to minimize this error, 
would be to split the line in more sections, as was attempted. 

Finally, we note that the current distribution does not show the 
marked peak near the resonance region obtained with parallel feed type of 
structures, or circuits. On one side this seems a reasonable result due 
to the different physical lay-out of the antenna; but, on the other side, 
it remains to be explained why the radiation patterns show acceptable 
results with different current distributions. An inspection of the phase 
distribution of the current (column GANG) shows that in all cases the 
current at the feed point of any two consecutive elements is very nearly 
in phase opposition, except for the element with dimensions closest to the 
resonance dimension, and the one before that. For these, the phase 
difference is near 90°, as also happens in other LP structures. This means 
then that only the **active region** is important for radiation characteristics, 
and the two models of current distribution previously used are essentially 
equivalent to each other and to the computed values. 
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APPENDIX A 



!• Geometric Relations in the BLPZZ 

With the aid of Figo A-1 the following relations are evident: 



L = — Y 
n n-1 



and 



b = 

n n-1 



(A-l) 



The relation between the length of the radiating wires and the stubs 
can be obtained in the following way: 



L = 



2b^ cosy 



2b’ cosy 



I _ 



n 



n tan aE 



n tan as 



(A-2) 



2b cosy 2b - cosy 

L ~L = — r = 2b ^siny + 2b siny + 4b ^ ’siny 

n n-1 ' tan aE tan aE n-1 ' n n-1 ' 



cosy b ^ 1 

— ^ — ( — - 1) = b - siny (1 + — ) + 2b / siny 

tan aE p n-1 ’ p n-1 



(l-p) = tanr (1+p) + 2b^_^' p tany 



_n 2p tany tan aE 



b^’ (l-p) - (1+p) tany tan aE 



(A-3) 



Other useful relations are 



2L’ - = L 1 + L ^ tan aE tany + L -L tan aE tany 
n-1 n-1 n-1 n n 



L .(!+— )-hL -(1- — ) tan aE tany 
n-1 p n-1 p 



= n-1 



[(1+p) + (p-1) tan aE tany] 
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^ L 

f 




a) Basis structure 



V 




b) Bent structure 

Fig A-1 - Geometric Relations in the BLPZZ 
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from which 



b' = [(1+p) - (1-p) tan aE tanv] tan as 

n-1 4p COSY 



and combining with (A-3) 

_ 2/3 tan‘^£. 

4 tacosy 



and, since 



n-1 



2b , cosy 

n-1 



tan aE 



(i-fi) ~(h-Id) tan^ tano^£ = {Up) iar)^ta'no(s-{l-p‘)t^'>^'^£.tar)0!siar)^^ 






^ J-tAj I" tanc>(£ + tancK^) - J (~r^) tar?oi'£ tarycx^ 



2 L 



»)0(E 



Ur>o(^ 



(A-i) 



Since the structure is usually defined by the parameters aE, as, i, 
b for the last element and the number of elements, equation (A-4) gives 
the next parameter, y> equations (A-1) and (A-2) define the rest of 

the structure. The length of the stubs can be obtained with (A-3), or 
approximately by 



b* 



n 



n tan as 
r' tan aE 

/p 



(A-5) 



This equation is exact when as = aE, but it can be used as a good 
approximation for the typical values used. The equation is obtained by 
inspection of Fig. A-1, with the help of a sort of **optical illusion." 

In fact, when as = aE the structure repeats itself at increments, with 
an upside down rotation. When as ^ aE, however, the rate of change of b’ 
with as is not linear and the equation is not true. 

The relations concerning the "phase center" of the element are 
obtained from Fig. 3, by inspection 
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(A-6) 



(oTl siyy £ i- cosCi^~S) 

6 >, Sin (}f-t) 



lav ^2 - si'n E ■*■ ^ 



in (r^) 



(A- 7) 



^1 = 



in smia-l) 
Cos S) 



(A-8) 



Cj - ^r? 



Cos 






(A-9) 



II . Transformations of Coordinates 

Expressions for the angle between a wire and a direction -6’. 
Considering Fig* A-2, that is a planif ication of the Y=0 and Z=0 planes, 
we obtain the following relations: 

cos a sin b = cos0 

cos a cos b = sin6 cos 0^^ 0^^ = 180-0 

R COS0’ = R cos a cos (0-b) 

= R [cos 6 cos a cos b + sin 6 cos a sin b] 

= R [cos B sin 0 cos 0^^ + sin B cos 0] 

= R [sin B cos 0 - cos B sin 0 cos 0] 

The wire in the figure corresponds to wire 1 in a BLPZZ element, and 
the angle B is given by (see Fig. A-1) 

B = e + 90® - Y = 90° - (y-e) 
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Fig, A-2 - Transformation of coordinates 
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so that, for wire 1 



Cos $1 ^ cos {^-i) C0S$ -Sir) sir) $ cos ^ 



(A-10) 



For wire 2 we have instead of an angle 65 and angle -65 given by 



-6 = - (90 - (Y+e)) 



and using the equation for cos we get 

cos ^ Cos cos e - sin (gie) Sm 0 cos ^ 



(A-11) 



For wire 3, we obtain (see Fig. A-3) 
cos a sin b = cos 0 
cos a cos b = sin 0 cos 0 

« 

cos 0 ^ = cos a cos ( 6 ~b) 

= cos a cos b cos 6 + cos a sin b sin B 

= cos 0 sin B + sin 0 cos B cos 0^ 

and, with 0^ = 0 and B = 90 - (y-e) 

cos0^ = cos 0 cos (y"c) + sin 0 sin (yc) cos 0 (A-12) 

For wire 4, we note that the only difference is again the direction 
of the angle B = 90 - (y+s)> so that 

cos0^ = -cos 0 cos (y+s) + sin 0 sin (y+c) cos 0 (A-13) 
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Fig. A-3 - Transformation of Coordinates - Wire 3 
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III. The Components of the Vector In Spherical Coordinates 

According to Fig. A-4, we obtain the following relations, for wire 1 



^0, -- ^e; 


cos /4, 


(A-14a) 




sin ,A| 


(A-14b) 



Using the expression from spherical trigonometry 



cos A^ 



cos a - cos b cos c 
sin b sin e 



we obtain 



cosA, - cos(^o-A)-cosd,'cos$ _ imp- cose (sir)pcos$-cosps/n0cos^) 
Sm 6) 5/0 6 \J I- cos^ 0^' 5/Vj S 



smp~ s'lnpcos^e f c.osp sinB cos$ cos p s,n pO~cos^$) +cosps,inOcosdcosp 

^ l-cos^; Jhcos^ej 



5/r7 P Sin $ f cosp Cos e COSp 

)/ h cos^O\ 



and, since g = 90 - (y-e), 

cos A,. CoS Of-e) Sir) B * sir) (ps) cos6co$ ^ 

\j /- cos^o] 



(A-15a) 



We have the denomination in the form 9, because in the field 

expression we have a sin 0’ in numerator that cancels with this. Now we 
obtain, from here, sin A 
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Fig. A-4 - Components of the field vector 
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Sin^A, = I-COS^A, - !~( ^os^/B Sin^d co<.^fl!-2sirtficos^ $m$cos6 cospf) 



sm^9 * CO^A coi^^tz $m[icosp SM$Cosdcos^) 
SfY) ^ 6 \ 



that reduces to 



^ J-Siy^^p-Ccs^f/SO'Sin^) ^ (h sin^/^) (hCcS^^) ^ cos^/8 S/r?^fi 

s'm^Ol * sm^0l ”* sm^Ol 

and 



S/r^ ^ COSjS ^ 

^>/V) ($?'' 



Sirt(^C-i) 5/V) ^ 
Sir, (>l 



The procedure is analogous for the other wires, and the results 
for wire 2 






■ Eg. sin Aj 



COS A 



2 z 



CoS Ori-£) sir) 9- Sm(Kte) cos9 cos\ 
sir, 0 / 



Sir) Aj - •5'’^ 



Sin $[ 



(A-15b) 

are 

(A-16a) 

(A-16b) 

(A-17a) 

(A-17b) 
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for wire 3 



E 



0 



3 




cos 



(A-X8a) 




sin A3 



(A-18b) 



CosA^ _ Cos (g’-e ) sin Siy^ Of- e) cos Ocos^ 



Sir) 



(A-19a) 



SinA^ - 



sir)(y-e) sin gf 

Sir,0^ 



(A-19b) 



and for wire 4 



\ cos A^ 



(A-20a) 



\ = ^e; ^ 

4 4 



(A-20b) 



cos A^ ^ cos(g+e) sine <■ sin (g+i) cos9 cos ^ 
sin 



(A-21a) 



sin _ s/W<'r^ej Sm (A-21b) 

Sin Oc, 

IV. The Angle a for the Correction of Reference Point 

This is the angle between the direction of c (see Fig. 3) and the 
direction to a point in the far field. Since c is in the same vertical 
plane as the wires, we can use the relations developed for the conversion 
of the angle 0* to spherical coordinates, with 6 in place of 3, and we 
obtain the following relations: 



cos 


“l = 


- [sin 




COS 


0 


+ 


cos 




sin 


0 


cos 


0] 


(A-22) 


cos 


“2 = 


sin 




cos 


0 


- 


cos 




sin 


0 


cos 


0 


(A-23) 


cos 


“3 = 


sin 


1— 1 

CO 


cos 


0 


- 


cos 




sin 


0 


cos 


0 


(A-24) 


cos 


II 


sin 


CM 

CO 


cos 


0 


4- 


cos 


^2 


sin 


0 


cos 


0 


(A-25) 
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,VKC 



/<? 







iviVe i 














c) v^*«<e ^ 



C 




^ ) v^)Ve 4 



Fig. A-5 - Reference point and direction angles 
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The angles 
(A- 6 ) and (A-7). 



and 62 are shown in Fig. 3, and given by equations 
A-5 shows the way how the a angles are defined. 



68 



APPENDIX B 



Field Equations for Wires 2, 3 and 4 
I - Wire 2 

Since wire 2 (see Fig. 5) is in series with wire 1, we need to add in 
the exponent a factor corresponding to the length of wire 1 - 2b - and we 
get from (7) 



r. e 






e - e 






(B-l) 



Now, using equations (A-11), (A-16) , (A-17) and (A-23) , we obtain 



I. e 



(B-2) 






‘ ^ ^ 

[o< + ^ k - C<>5 j 

^COS(gre) Sir) & ~ S/n (^f£) COS $ cos/j 

-jk[3-|bfC2C0^of^] [o0^k(^~cos6>J)]L 



Sirt (^*-e) s/y> / 
correction foi 

being C 2 and cos U 2 given by equations (A-9) and (A-23). 



— ik c cos ci 

where the factor e 2 2 is the correction for the reference point, 



II - Wire 3 

Here we consider the current flowing upward, but to use the same 
exponential representation for the current, we consider the distance as 
negative away from the source. This is the usual representation for trans- 
mission lines, inherent in the present assumption of a travelling wave in 
the antenna. Then 
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and 



^3 ^0 ® 



-jBs 



E., ^ ^ 1, 5,. 0' e e 






j:_£_ 






(B-4) 



Now, using equations (A-12), (A-18), (A-19) and (A-24), we obtain 






[ |. 6 f C3 coscK^J ^ij] b - [o< +jl< ^i)] b 

X [^<’5 5/V) (9 - S/n (/'£) cos 9cos^] (B 5) 



_^k[-|-b+^3^®-S‘=^i] -[(^+^/<('|:K^s^i)Ji> 

hs25^=-HI, e _£ ^ :,_£ 



X 5/n (^-e) sin JZ^ 



X 

(B-6) 



where c^ and a^ are given by equations (A- 8 ) and (A-24). Note that 
C 3 = and (Fig. 5). 

III. Wire 4 

This wire is in series with wire 3, so that the procedure is analogous 
to wire 2 that was in series with wire 1 

.jpb [o(+jl<(^fC05 6i.'j]i5 _[o<+jlc(-|-+t«:S<?4)]b 

E. ^ £k L 5m 0j £ :_g (B-7) 



Now, using equations (A-13) , (A-20) , (A-21) and (A-25), we obtain 
.Lf^3b+C4C05c<^,'] [„( +^I<(t|. +C 0iZ?^'j]b -[c<*i^i-^+<^os&^)]b 

I ^ ^ >= 






X Si'rid + 5m Q^■d) cosd cos^j 



(B-8) 
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-fx+jK-f- 

> --_£ 



X sin Q^£) .5(V) 



&J,)]L 

(B-9) 
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APPENDIX C 



Equations for Non-Uniform Line 
I. Line Parameters 

Neglecting wire resistance and conductance between wires, the line 
parameters are 



Z = j 120 k In (s/a) 



and 



Y = j k/ (120 In (s/a)) 
as given by equations. 

From here we obtain the line impedance 



K = /zTy = 120 In (s/a) = -fr 

J ^ 



and the propagation factor not accounting for losses 



r = /zY = jk 



In terms of the structure defining dimensions and angles (see Fig. 
A-la) , we obtain the following equations 
For the radiators 



Z(x) = j 120k (B + In (A + x)) 
Z(x) = j 120k (C + In (D - x)) 



for 0 _< X < 2b (C-la) 
for 2b < X < 4b (C-lb) 



where 



1 - 2b siny 
cos (y-e) 



sine 



B = In 



2 cos (y-e) 
a 
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C = In 



2 cos (y-He) 
a 



L " 2b sin y 
cos (Y+e) 



sinc+ 



2b cos (Y-e) ^ 2b 
cos (Y+€) 



where 



For the stubs 

z(x) = jl20K(^ E + In (F+oc) ^ 



(C-2) 



_ , 2sin vsin e 

E = In ‘ 



2b tan a 



F = 



sin Y JF a 



and the quantities involved are as defined in Fig. A-la 




Fig. C-1 Input Impedance. 
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!!• Radiation Resistance 

According to Chaney and Klamm [24], [25], the radiation resistance of 
rhombic antenna in free space is given by 

R^ = 120 1^1.1544 + 2 In (2 kl sin^ a) + 2 Ci (2kl) - 2 Ci(2kl sin a) 

- Ci (2kl (1 + cos a) - Ci (^2kl (1 - cos a) ^ 

+ cos ^ 2kl sin^a ^ j^Ci (^2kl cos a (1 - cos a) ^ + Ci (^2kl sin a (1 - 

+ Ci (^2kl cos a(l+cos a)^ + Ci^2 kl sin a(l+sin a)^ - 2Ci ^ 2 kl cos^ a ^ 



- 2Ci Q 2kl sin^ a ^ 



;in (\ 



+ sin ( 2kl sin a 



0 [ ^ (1-cos a)^ - Si (^2kl sin a (1-sin a)^ 



- Si ^2kl cos a (1+cos a)^ + Si ^2kl sin a (1+sin a) ^ + 2 Si 2kl cos^ 



< 



- 2 Si ( 2kl sin a 



where 1 is the wire length - 2b - , a is the vertex angle 90 - Y > 

Ci is the cosine integral 



ci(x) = J 



COS V 



dv 



sin a) 



) 

)]} 
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and Si is the sine integral 



Si(x) 



/ 



0 



sin V 

V 



dv 



III . Input Impedance and Current Distribution 

Using equations (33) and dividing the line in sections as shown in 
Fig. D-1, we obtain for each section 



Z. 

in 



(I,n) 



KaveCi n') Kave(I.n) B(I,n) + Zin(I.n+l) D(I,n) 
Kaveu.n; 2in(I,n+l) C(I,n) + Kave(I,n) A(I,n) 



Each stub counts as only one section of uniform line. The number M 
of divisions within each radiation depends on the degree of approximation 
desired, but it is convenient to use for M an even number, to avoid the 
relatively faster rate of change of characteristics at the middle point. 
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APPENDIX D 



Computer Programming 
I. Program AZORES 

This program is written in code FORTRAN IV and it was prepared for 
the IBM 360 system that is operating at the Naval Postgraduate School. 

The program uses two subroutines. The main program reads the data, 
computes and writes the physical characteristics of the antenna, selects 
the current distribution for the antenna wires, controls the coordinates 
of the desired pattern, prints the results, and calls for the subroutines. 

The pattern computation is made by subroutine CORVO, that is given 
the antenna parameters, the elevation and azimuth angles, and the current 
distribution. In output, CORVO gives the two components of the electric 
field. 

The program computes each pattern twice, for two different models of 
current distribution. Model 1 (see Fig. 9) is fed in data cards; Model 2 
is generated by the main program. 

The use of a large number of comment cards, hopefully makes the 
program self-explanatory. As much as possible, the names given to 
algebraic and geometric symbols have the same phonetics as the ones used 
in the text. The only remarkable exception is the s 3 rmbol TAU, that stands 
for p, the structure scaling factor. 

The following pages contain the programming statements and the output 
print of a sample program. 
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SAMPLE PRCGRAM AZORES - MAIN 
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SAMPLt PRCGPAP AZCRES - MAIN 
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ALFA2 =-(SIN(OELTA2)^COS(TETA)-COS(DELTA2)^SIN(TETA)^COS(Fn ) 
ALFA2A= ( S IN (CELT A2) ♦COS (TETA) +COS( DELTA2) ♦SIN ( T ETA ) ♦COS ( F I ) ) 
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-90.000 


-13. 591 


C.O 


65. 


0. C8A 


-90.000 


-18. 133 


0 . 0 


70. 


0.033 


-90.000 


-26. 299 


0 . 0 


75. 


C.Cll 


90.030 


-36.072 


!'• I 


CC. 


C . CA6 


9C.000 


-23. 394 


n 


85 . 


C . G 7 3 


90. 000 


-19.352 


0 . c 


90. 


C.09 3 


9C. 030 


-17.244 




C5. 


0. 1C 7 


90.000 


-16.003 


0.0 


IOC. 


C. 117 


90.000 


-15.221 


C .0 


105. 


C. 12A 


9C.0Q0 


-14.694 


0 . r 


1 1C. 


C. 130 


90.000 


-14. 306 


P 

U.' • L. 


115. 


C. 13A 


90.000 


-13. 997 


0.^ 


12C. 


C. 139 


90.000 


- 13. 741 


0 . 0 


125. 


0. 1A2 


90.000 


-13. 535 


0.0 


130. 


0. lAA 


90.000 


-13.334 


. c 


135. 


C. 1A6 


90.000 


-13.294 


0 . 0 


1 AC. 


C . 1A6 


60. 000 


-13. 258 


c.c 


1A5. 


C. 1A6 


9C. 000 


-13. 301 


0.3 


150. 


C. lAA 


9C. OCO 


-13.381 


0.0 


155. 


C. 1A3 


90.000 


-13.493 


C.C 


16C. 


C. lAO 


9C. 000 


-13.618 


0 . 0 


165. 


0.139 


9C.C00 


-13.739 


C.O 


170. 


C . 1 3 7 


90.000 


-13.839 


c.c 


175. 


C. 136 


90.000 


-13.904 


0 . sj 


IRC. 


C. 136 


90. COO 


-13.927 


UIRFCTIVITY/ISC 


FI 3DB 


ELEV 30B 




1 0 . 9 6 9 c e 


3^. 


25. 







THIS CIF.F.CT1VITY FIGURE V^AS COMPUTED 

ASSUMING THAT THF PATTERN IS VvEL L-8EHAVFD . 
IT HAS NO MEANING IF THERE ARE SIDE-LOBES. 
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CURRENT 


niSTRIBUTICN 


2 


FREO 


= 0.50CF 09 




ELEV 


F I 




E 


PSI 


EDR 


0.0 


0.0 




0.660 


-90.000 


- 0 . 0 oc 


10. 


C.C 




0.6A? 


-9C.OOO 


-0.348 


2C . 


0.0 




0.56 7 


-90.000 


-1.433 


30 . 


0.0 




C.A52 


-90.000 


-3.393 


AC. 


c.c 




C. 315 


-90. COO 


-6. 529 


50 . 


c.c 




0. 179 


-90,000 


-11.434 


60 . 


o.c 




0.071 


-90.000 


-19.472 


70 . 


C.C 




C.CIO 


-90.000 


-36.778 


ec . 


o.c 




0.CC7 


90.000 


-39.533 


<5C . 


leo. 




C. 002 


-9C.000 


-52. 310 


80. 


180. 




C.CCA 


90.000 


-43. 891 


70. 


IfiO. 




C.C09 


90.000 


-37. 003 


60. 


ISC. 




C.CPfl 


90.000 


-27. 481 


5C. 


leo. 




0.069 


90.000 


-19.740 


AC. 


180. 




0.122 


90.000 


-14.736 


30. 


180. 




0.172 


90.000 


-11.764 


20. 


180. 




C. 208 


9C.OOO 


-10.153 


IC. 


1 80. 




C.226 


90.000 


-9.403 


0 .0 


180. 




C.232 


90.000 


-9. 191 


■N ^ 


C.C 




C.668 


-90.000 


0.0 


C .0 


5. 




C. 663 


-90.000 


-0.062 


C .0 


10. 




r . 6A9 


-90.000 


-0.249 


0.0 


15. 




0.625 


-90.000 


-0.573 


r n 


20. 




C. 59 2 


-90,000 


-1.050 


0.0 


25. 




C.5A9 


-90.000 


-1.705 


0.0 


30. 




0.AO7 


-90.000 


-2.578 


0.0 


35. 




C.A35 


-90. 000 


-3. 724 


o.c 


AO. 




C .366 


-90.000 


-5.226 


c.c 


A5. 




C. 29 1 


-90.000 


-7. 213 


0 . 1,. 


50. 




0.213 


-9C.000 


-9.914 


c.o 


55. 




0.136 


-90.000 


-13.815 


C . 0 


60. 




C.C6A 


-90.000 


-20.413 


c.c 


65. 




G.CCO 


-89. 95A 


-75.435 


0 . 0 


70. 




C.C51 


90.000 


-22.317 


0.0 


75. 




0. C88 


90.000 


-17.601 


o.c 


80. 




0.110 


9C.000 


-15,636 


u . C 


85. 




0.120 


9C.000 


-14, 907 


0.0 


90. 




C. 121 


90.000 


-14. 863 


'I ry 

• \ i 


95. 




0.117 


9C.000 


-15.145 


r ** 

O . ’-r 


ICO. 




0.113 


90.000 


-15.427 


0.0 


1C5. 




0.113 


90.000 


-15.424 


G . C 


no. 




0.119 


9C .000 


- 14. 985 


0.0 


115. 




0.131 


90.000 


-14.153 


o.c 


120. 




0. 1A8 


90.000 


-13. 103 


V-* ♦ 0 


125. 




C . 1 6 7 


9C.000 


-12.029 


o.c 


130. 




0.187 


90.000 


-11.066 


0.0 


135. 




C.20A 


90.000 


-10.285 


c.o 


lAO. 




0.219 


90.000 


-9.708 


c.o 


1A5. 




0.228 


9C.000 


-9.323 


0 * L 


150. 




0.23A 


90.000 


-9. 104 


c.o 


155. 




0.237 


90.000 


-9.013 


c.o 


160. 




0.237 


90.000 


-9.011 


O.Q 


165. 




0.235 


90.000 


-9.060 


0.0 


170. 




C.23A 


90.000 


-9. 12? 


c.c 


175. 




0.232 


90.000 


-9. 172 


c.o 


180. 




C.232 


90.000 


-9. 191 


DIRECT IV I TY/ I sc 


E 1 


3CB 


ELEV 3DB 





lC.969Ce 33. 25. 

THIS DIRECTIVITY FIGURE WAS COMPUTED 

ASSUMING THAT THE PATTERN IS WELL-REHAVED. 
IT HAS NG MEANING IE THERE ARE SIDE-LOBES. 
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CVFkALL LtNGTF CF STHLCTURE 3.541 



CURRENT D 


TSIRinUTICN 1 


FRFO= 


G.llOF 10 




ELFV 


F I 


F 


PSI 


FDB 


O.C 


c.c 


3.228 


9C.OOO 


-0.000 


1C . 


C.O 


3. 1A3 


90.000 


-0.232 


20. 


0 . 0 


2.9C0 


90.000 


-0.93C 


3C . 


c.c 


2.537 


90.000 


-2.093 


AO . 


c.o 


2.10A 


90.000 


-3. 7?0 


5G . 


C .0 


1.653 


00,000 


-5.816 


6C . 


C.( 


1.225 


90. OOO 


-8. A16 


7C . 


0.0 


0.8A6 


90.000 


-11.636 


80 . 


c.o 


0.525 


90.000 


-15. 768 


90 . 


180. 


C.270 


-9C.000 


-21.539 


80. 


180. 


C.C82 


-89.999 


-31.878 


70 . 


IPC. 


C.0A3 


89.999 


-37. A85 


60. 


180. 


0.117 


90.000 


-28. 789 


50 . 


lec. 


C. 156 


90.000 


-26.308 


AO. 


IHO. 


C.17A 


•RC.OOC 


-25.358 


30. 


1 80 . 


0 .182 


90.000 


-2A.97A 


20. 


1 PC. 


C . 186 


90.000 


-24. 793 


10. 


180. 


C. 1 88 


9f .000 


-24.690 


c .c 


IHC . 


C. 189 


9C. OCO 


-24.653 


C . 0 


C.c 


3.228 


9C.000 


0.0 


O.C 


5 . 


3.22A 


9C.000 


-0.012 


0.0 


10. 


3.211 


90.000 


-0.047 


0.0 


1 5 . 


3.189 


90.000 


-0. 106 


0 . 0 


20. 


3.158 


9C.00O 


-0. 190 


0.0 


25. 


3.119 


9C. 300 


-C. 300 


0.0 


3C. 


3. 070 


9C .000 


-0.437 


O.C 


35. 


3.311 


90.000 


-0.605 


C . 0 


AC. 


2.9A 1 


90.000 


-C. 808 


,'■> 

w . o 


A5. 


2.860 


90.000 


-1.05? 


0.0 


5C. 


2.765 


9C.030 


-1.344 


O.C 


55. 


2.656 


90.000 


-1.694 




6c. 


2.531 


90.000 


-2.114 


c.o 


65. 


2.389 


9C.00C 


-2.616 


c.o 


70. 


2.230 


90.0CG 


-3.213 


O.C 


75. 


2.056 


9C.900 


-3. 920 


c.o 


80. 


1 . 869 


9C.00C 


-4. 748 


0.0 


S5. 


1.673 


9C.C03 


-5.707 


O.C 


9C. 


1 .A75 


90.000 


-6.802 


0.0 


95. 


1.280 


9C.000 


-8.034 


O.C 


100. 


1.C9A 


90.000 


-9. 397 


O.C 


1C5. 


C.923 


90.000 


-10.877 


c.o 


110. 


C .77C 


90.000 


-12.449 


0 . 0 


115. 


0.63 8 


9C.000 


-14.079 


c . c 


120. 


C.528 


9C.000 


-15.722 


c.c 


125. 


C.A39 


90. COO 


-17.325 


c .0 


130. 


C. 369 


90.000 


-18.829 


O.C 


135. 


C.316 


90. OCO 


-20. 184 


0 .0 


lAO. 


0.276 


90.000 


-21 . 350 


0.':' 


1A5 . 


C.2A7 


9C.000 


-22.310 


c.c 


15C. 


0.227 


90.000 


-23.064 


O.C 


155. 


0.212 


9C.000 


-23.633 


0.0 


160. 


C.2C3 


90. 000 


-24.046 


c.o 


165. 


C. 196 


90.000 


-24.331 


0.0 


170. 


C. 192 


90.000 


-24.516 


c.o 


175. 


C.19C 


90.000 


-24.619 


O.C 


180. 


0. 189 


9C.000 


-24.653 


GIRECTIVITY/ISO FI 


3DB 


FLEV 308 




6.368DG 


68. 


35. 







IHIS DIRECTIVITY FIGURE UAS Cn^^PUTE^ 

ASSUMING THAT THE PATTERN IS VvEL L- BEHAVED . 
IT HAS NO i^EANING IF THERE ARE SIDE-LOBES. 
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CURRLNT r IS ru IBUT ICN 2 FPEU= O.llOE 10 



El. LV 


F I 


E 


PSI 


EOB 


G .0 


O.C 


2.951 


90.000 


0.0 


1 


o.c 


2.87Q 


90.000 


-0.216 


20. 


c.o 


2.675 


90.000 


-0.852 


30. 


f'.c 


2.377 


90.000 


-1.878 


^0 . 


o.c 


2.C29 


90.000 


-3.254 


50. 


c.o 


1.664 


9C.000 


-4.976 


6C . 


c.c 


1.297 


9C.000 


-7. 14? 


7C. 


c.o 


C.93C 


9C.C00 


- 10.0 28 


. 


c.o 


0.580 


9C.000 


-14. 1 32 


9C . 


180. 


C. 284 


-90.000 


-20. 345 


80. 


18C. 


0.079 


-89.999 


-31.488 


70. 


180. 


C.C31 


89.999 


-39.676 


60. 


180. 


C.C78 


90.000 


-31.517 


50. 


18 0. 


C. 1C9 


90.000 


-28.629 


4C . 


180. 


C. 149 


90.000 


-25.949 


30. 


180. 


C. 197 


90.000 


-23.494 


20. 


180. 


C.243 


9C.OOO 


-21.678 


1 0 . 


180. 


C ,215 


90.000 


-20.611 


C .C 


180. 


C.286 


90.000 


-20.264 


0 . 0 


C.O 


2.951 


90.000 


0.0 


C .0 


5. 


2.949 


90. 0'''0 


-0.096 


C .0 


10. 


2.943 


9C. 000 


-0.024 


0 .0 


15. 


2.934 


90.000 


-0.052 


C .0 


20. 


2.921 


90.000 


-0.088 


0.0 


25. 


2.906 


90.000 


-0.133 


G . G 


30. 


2.889 


90.000 


-0. 184 


0.0 


35. 


2.869 


9C.OOO 


-0.244 


0.0 


40. 


2.845 


9C.O00 


-0. 317 


0 . 0 


45. 


2.814 


90.000 


-0.413 


c 


50. 


2.772 


90.000 


-0.544 


C . 0 


55. 


2.713 


9C.000 


-0.731 


c . c 


60. 


2.632 


90.000 


-0.996 


c .c 


65. 


2.522 


90.000 


-1.365 


0 .0 


70. 


2.380 


9C.OOO 


-1.867 


c.o 


75. 


2.205 


90.000 


-2.531 


c.c 


80. 


1.999 


90.000 


-3.384 


0 .0 


85. 


1.768 


9C.OOO 


-4.448 


o.c 


90. 


1.524 


90.000 


-5.742 


c.o 


95. 


1.277 


90.000 


-7.275 




IOC. 


1.04 2 


90.000 


-9.045 


c.c 


105. 


0.829 


90.000 


-11.030 


0 . 0 


110. 


0.64 7 


90.000 


-13.177 


0.0 


115. 


0.5C2 


90.000 


-15. 388 


o.c 


120. 


C.393 


90.000 


-17.512 


c.c 


125. 


C.318 


90.00C 


-19.351 


o.c 


130. 


0.272 


90.000 


-20.720 


'o • 


135. 


0.246 


9C.000 


-21.525 


0 . 0 


14 0. 


0.239 


90.000 


-21.817 


c.c 


145. 


C. 241 


9C.OOO 


-21.749 


0.0 


15C . 


G. 24 9 


9C.00C 


-21.488 


o.c 


155. 


0.258 


90.000 


-21. 164 


c.o 


16C. 


0.26 7 


90.000 


-20.853 


0.0 


165. 


C.275 


9C.000 


-20. 597 


. c 


17C. 


0.281 


90.000 


-20.412 


c.o 


175. 


C. 285 


90.000 


-20.301 


■> p 

-V • J 


18C. 


C.286 


9C.000 


-20.264 


LIRECT IV I TY/ ISO 


FI 3DB 


ELEV 308 




5. 772CP 


78. 


35. 







THIS DIPECTIVITY FIGURE UAS COMPUTED 

ASSUMING THAT THE PATTERN IS UELL-BEHAVED, 
IT HAS \C 'LEANING IF THERE ARE SIDE-LORES. 
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OVERALL LENGTH OF STRICTURE 0.541 



CURRENT 


CISfRlPUTICN 1 


FREQ = 


0.170E 10 




ELEV 


F I 


E 


PSI 


EOB 


O.C 


O.C 


2.330 


90.000 


0.0 


1C . 


c.o 


2.257 


90.000 


-0.275 


20. 


O.C 


2.052 


90.000 


-1. 103 


30. 


C . 0 


1.749 


90.000 


-2. 493 


40 . 


0.0 


1.394 


90.000 


-4.462 


50. 


O.C 


1.035 


90.000 


-7.051 


60. 


O.C 


0.708 


90.000 


-10. 344 


70. 


0.0 


C.437 


90.000 


-14. 533 


80 . 


O.C 


0.233 


90.000 


-20.013 


90. 


180. 


0.C97 


-90.000 


-27.627 


BC. 


180. 


C.025 


-90.000 


-39.335 


70. 


180. 


0.005 


-89.999 


-53.477 


60 . 


180. 


0 . C 1 6 


-9C.000 


-43. 334 


5C . 


180 . 


0.035 


-90.000 


-36.506 


40. 


180. 


0.045 


-90.000 


-34.311 


30 . 


180. 


0.040 


-90.000 


-35.350 


20. 


180. 


0.024 


-90.000 


-39.586 


10. 


180. 


C.C0 9 


-9C .000 


-48.263 


0.0 


180. 


C.CC3 


-90.000 


-58. 853 


O.C 


C.O 


2.330 


90.000 


0.0 


0.0 


5. 


2.322 


90.000 


-0.029 


O.C 


10. 


2.299 


90.000 


-0.118 


O.C 


15. 


2.259 


90.000 


-0.269 


O.C 


20. 


2.203 


90.000 


-0.486 


O.C 


25. 


2.132 


90.000 


-0.773 


0.0 


30. 


2.043 


90.000 


-1.140 


O.C 


35. 


1.939 


90.000 


-1.595 


O.C 


4C. 


1.819 


90.000 


-2. 149 


0.0 


45. 


1.685 


90.000 


-2.817 


0 . c 


50. 


1.537 


90.000 


-3.614 


0 . 


55. 


1.378 


90.000 


-4.560 


C .0 


60. 


1.212 


90.000 


-5. 630 


0.0 


65. 


1.040 


90.000 


-7.005 


0 . 0 


70. 


C.868 


90.000 


-8.575 


0.0 


75. 


0.700 


90.000 


-10.445 


O.C 


80. 


0.540 


90.000 


-12.699 


c.o 


85. 


0.392 


90.000 


-15.474 


Q.O 


90. 


C. 261 


90.000 


-19.029 


O.C 


95. 


0.147 


90.000 


-23.981 


O.C 


ICO. 


0.054 


90.000 


-32.655 


O.C 


105. 


0.018 


-9C.001 


-42. 158 


O.C 


1 1C. 


0.071 


-90.000 


-30.357 


0 . G 


115. 


0.105 


-90.0Q0 


-26.922 


0.0 


120. 


0.12 3 


-90.000 


-25.516 


0.0 


125. 


0. 129 


-90.000 


-25.138 


0.0 


130. 


0.125 


-90.000 


-25.442 


O.C 


135. 


0.113 


-90.000 


-26.278 


0.0 


140. 


0.097 


-90.000 


-27.579 


O.C 


145. 


0.080 


-90.000 


-29. 327 


0 . c 


1 50. 


0.062 


-90.000 


-31.547 


O.C 


155. 


0.04 5 


-90.000 


-34.306 


0.0 


16C. 


0.030 


-90.000 


-37.736 


c.o 


165. 


C . C 1 8 


-90.000 


-42.058 


O.C 


170. 


C.CIO 


-90.000 


-47.606 


O.C 


175. 


0.004 


-90.000 


-54.423 


O.C 


1 «0. 


0.003 


-90.000 


-58.852 



DIRECT IVITY/ISC FI BOB ELEV BOB 

7.881CR 48. B5. 



THIS DIRECTIVITY FIGURE WAS COMPUTFD 

ASSU'^INC THAT THE PATTERN IS WELL-BEHAVED. 
IT HAS NO NEANINC IF THERE ARE SIDE-LOBES. 



95 



1 01 ^ 



CURRENT niSTRIBUTICN 


2 FREQ= 


0.170E 10 




ELEV 


F I 


E 


PSI 


EDB 


O.C 


0.0 


2.235 


90.000 


0.0 


10. 


0.0 


2.176 


90.000 


-0.233 


20. 


O.C 


2.009 


90.000 


-0.927 


30. 


0.0 


1.760 


90.000 


-2.079 


^0. 


c.o 


1.A61 


90.000 


-3.693 


50. 


0.0 


1. lAA 


90.000 


-5.820 


60 . 


0.0 


0.831 


90.000 


-8. 598 


70. 


0.0 


0.5A1 


90.000 


-12.327 


80. 


O.C 


0.295 


90.000 


-17.600 


90. 


160. 


C. 116 


-90.000 


-25.686 


80. 


180. 


C.023 


-90. 000 


-39,767 


70. 


180. 


O.CIA 


-90.000 


-43. 819 


60. 


180. 


0.C66 


-90.000 


-30.559 


50. 


180. 


0.140 


-90.000 


-24.072 


AO. 


180. 


C.2CC 


-90.000 


-20.971 


30. 


180. 


0.229 


-90.000 


-19. 789 


20. 


180. 


0.231 


-90.000 


-19.724 


IC. 


18C. 


0.221 


-90.000 


-20.101 


0.0 


180. 


0.216 


-90.000 


-20. 317 


0.0 


0.0 


2.235 


9C.000 


0.0 


0.0 


5 . 


2.231 


90.000 


-0.019 


0.0 


10. 


2.216 


90.000 


-0.075 


0.0 


15. 


2.192 


90.000 


-0. 171 


0.0 


20. 


2.157 


90.000 


-0. 311 


0.0 


25. 


2.111 


90.000 


-0.499 


O.C 


30. 


2.052 


90.000 


-0. 742 


0.0 


35. 


1.981 


90.000 


-1.049 


0.0 


AO. 


1.895 


90.000 


-1.432 


0.0 


A5. 


1.795 


90.000 


-1.907 


0.0 


50. 


1.678 


90.000 


-2.490 


0.0 


55. 


1.5A5 


90.000 


-3.206 


0.0 


60. 


1.397 


90.000 


-4.081 


0.0 


65. 


1.235 


90.000 


-5.153 


0.0 


70. 


1.061 


90.000 


-6. 470 


0.0 


75. 


0.880 


90.000 


-8.099 


0.0 


80. 


0.695 


90.000 


- 10. 144 


0.0 


85. 


C.513 


90.000 


-12.787 


O.C 


90. 


0.338 


90.000 


-16.404 


0.0 


95. 


0.177 


90.000 


-22. 045 


c.o 


100. 


0.033 


90.000 


-36. 655 


0.0 


105. 


C.C9C 


-90.000 


-27.934 


0.0 


no. 


C. 189 


-90.000 


-21.466 


0.0 


115. 


0.26A 


-90.000 


-18.554 


0.0 


120. 


0.316 


-90.000 


-16.989 


0.0 


125. 


0.34 7 


-90.000 


-16. 173 


0.0 


130. 


0.360 


-90.000 


-15.850 


0.0 


135. 


C.359 


-90.000 


-15.880 


0.0 


lAC. 


0.347 


-90.000 


-16. 172 


0.0 


1A5. 


0.328 


-90.000 


-16. 657 


0.0 


150. 


0.306 


-90.000 


-17.276 


0.0 


155. 


0.282 


-90.000 


-17. 968 


0.0 


160. 


C.260 


-90.000 


-18.672 


0.0 


165. 


0.242 


-90.000 


-19. 322 


O.C 


170. 


0.227 


-90.000 


-19. 850 


0.0 


175. 


C.219 


-90.000 


-20. 196 


0.0 


180. 


0.216 


-90.000 


-20.317 


DIRECTIVITY 


/ISO 


FI 3DB 


ELEV 30B 




7. A510B 


53. 


35. 






THIS DIRECTIVITY FIGURE WAS COMPUTED 




ASSU*^ ING 


THAT THE 


PATTERN IS 


WELL-BEHAVED. 


IT HAS NO 


MEAN ING 


IF THERE ARE SIDE-LOBES. 
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ALL LMTS AKE <^KS 
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HEIGHT CF FEEC PCIM 0.003 




CURRENT 1 


D IS IR IRUT ICN 


1 


FREQ= 


0.200E 10 




ELFV 


F I 




E 


PSI 


EDB 


0 . C 


C.O 




2.581 


-90.000 


-0.000 


10. 


0.0 




2.A9A 


-90.000 


-0.297 


20. 


0.0 




2.2A9 


-90.000 


-1. 196 


30. 


0.0 




1.889 


-90.000 


-2.710 


AC. 


0.0 




1.A75 


-90.000 


-A. 861 


50. 


0.0 




1.067 


-OC.O'^O 


-7.672 


60. 


c.o 




C.71A 


-90.000 


-11.161 


70 . 


0.0 




C.AAO 


-90.000 


-15.363 


80. 


0.0 




0.2A7 


-90.000 


-20.38A 


9C. 


180. 




C. 120 


90.000 


-26.630 


80. 


180. 




C.OAO 


90.000 


-36, llA 


7C. 


180. 




0.012 


-90.001 


-A6. 7A0 


60. 


180. 




0.052 


-90.000 


-33.987 


50. 


180. 




0.089 


-90.000 


-29.259 


AO. 


180. 




0.12 8 


-90.000 


-26.067 


30. 


180. 




0.169 


-90.000 


-23.68A 


20. 


180. 




0.205 


-90.000 


-21.993 


10. 


180. 




0.231 


-90.000 


-20.982 


0.0 


180. 




C.2AC 


-90.000 


-20.6A5 


0.0 


0.0 




2.581 


-90.000 


0.0 


0.0 


5. 




2.572 


-90.000 


-0.032 


0.0 


10. 




2.5A3 


-90.000 


-0. 129 


0.0 


15. 




2.A96 


-90.000 


-0.293 


0.0 


20. 




2. A29 


-90.000 


-0. 527 


0.0 


25. 




2.3AA 


-90.000 


-0.836 


0.0 


30. 




2.2A2 


-90.000 


-1.225 


0.0 


35. 




2. 122 


-90.000 


-1.702 


0.0 


AO. 




1.987 


-90.000 


-2.275 


0.0 


A5. 




1.838 


-90.000 


-2.952 


0.0 


50. 




1.677 


-90.000 


-3. 7A3 


0.0 


55 . 




1.509 


-90.000 


-A. 660 


0.0 


60. 




1.337 


-90.000 


-5. 712 


0.0 


65. 




1.165 


-90.000 


-6.910 


0.0 


70. 




0.997 


-90.000 


-8.263 


0.0 


75. 




0.838 


-90.000 


-9. 776 


0.0 


80. 




C.691 


-90.000 


-11.A50 


0.0 


85. 




0.560 


-90.000 


-13.276 


0.0 


90. 




0.AA7 


-90.000 


-15.226 


0.0 


95. 




C.35A 


-90.000 


-17.2A8 


0.0 


ICO. 




0.281 


-90.000 


-19.2A9 


0.0 


105. 




C.228 


-90.000 


-21.092 


0.0 


110. 




0.191 


-90.000 


-22.608 


0.0 


115. 




0.170 


-90.000 


-23.6A9 


c.o 


120. 




0.160 


-90.000 


-24. 153 


0.0 


125. 




0.160 


-90.000 


-2A. 178 


0.0 


130. 




0.165 


-90.000 


-23.867 


0.0 


135. 




0.175 


-90.000 


-23.375 


0.0 


lAO. 




0.186 


-90.000 


-22. 826 


0.0 


1A5. 




0.198 


-90.000 


-22.298 


o.c 


150. 




0.2C9 


-90.000 


-21. 832 


0.0 


155. 




0.219 


-90.000 


-21. AA6 


o.c 


160. 




0.226 


-90.000 


-21. 1A3 


0.0 


165. 




0.232 


-90.000 


-20.917 


0.0 


170. 




0.236 


-90.000 


-20. 76A 


0.0 


175. 




0.239 


-90.000 


-20.67A 


0.0 


180. 




0. 2A0 


-90.000 


-20.6A5 


DIRECTIVITY/ISO 


FI 


3C8 


ELEV 3DB 





T.aaiop A8. 35. 

THIS DIRECTIVITY FIGURE V<AS COMPUTED 

ASSUMING THAT THE PATTERN IS V^ELL-BEHAVED. 
IT HAS NO MEANING IF THERE ARF SIDE-LOBES. 
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CURRENT 


CISTKIBUT K.N 


2 


FREQ 


= 0.200E K 




LLEV 


FI 




E 


PSI 


EDR 


C . C 


0.0 




2.2H8 


-QO.OOO 


-0.000 


10. 


0.0 




2.201 


-90.000 


-0.335 


2C . 


0.( 




1.959 


-90.000 


-1.351 


3G. 


c.{ 




1.607 


-90.000 


-3.071 


AO. 


0.0 




1.210 


-90.000 


-5.530 


50. 


0.0 




0.835 


-90.000 


-8.757 


60 . 






0.527 


-90.000 


-12. 7A5 


70 . 


o.c 




0.308 


-90.000 


-17. A09 


B ) . 


0.0 




0.171 


-90.000 


-22. 5A6 


90 . 


18C. 




C.C90 


90.000 


-28.079 


80. 


IRC. 




0.038 


90.000 


-35.573 


70. 


180. 




C.OlO 


-90.001 


-47. 134 


60. 


180. 




0.069 


-90.000 


-30.415 


50. 


180. 




0. 1A3 


-90.000 


-24.070 


AO. 


IRO. 




C.228 


-9C.OOO 


-20.013 


3 C . 


lec. 




0.31A 


-90.000 


-17. ?45 


20. 


IBO. 




0.388 


-90. 000 


-15.423 


10. 


1 80. 




C . A 3 7 


-90.000 


-14. 385 


C.P 


IHO. 




0. A5A 


-90.000 


-14.048 


G .0 


O.C 




2.288 


-90.000 


0.0 


0 . 0 


5. 




2.278 


-90.000 


-0.040 


0.0 


10. 




2.2A6 


-90.000 


-0.160 


0 .c 


15. 




2. I9A 


-9C.000 


-0.363 


0.0 


20. 




2.123 


-90.000 


-0.651 


0 . 0 


?5. 




2.032 


-90.000 


-1.031 


0 . c 


30. 




1.923 


-90.000 


-1.509 


c.o 


35. 




1 . 798 


-90.000 


-2.093 


0.0 


AO. 




1.659 


-90.000 


-2. 791 


o.c 


A5. 




1.509 


-90.000 


-3.615 


0.0 


50. 




1.351 


-90.000 


-4. 574 


o.c 


55. 




1 . 190 


-90.000 


-5.682 


0.0 


60. 




1.028 


-90.000 


-6.947 


O • w 


65. 




0.872 


-90.000 


-8.380 


0.0 


70. 




0.725 


-90.000 


-9.984 


c.o 


75. 




0.59? 


-90.000 


-11. 750 


o.c 


80. 




0.A75 


-90.000 


-13.651 


0.0 


85. 




0.379 


-90.000 


-15.625 


c.o 


90. 




0.303 


-90.000 


-17. 555 


c.o 


95. 




0.2A9 


-90.000 


-19. 257 


o.c 


100. 




0.216 


-90.000 


-20. 506 


c.o 


105. 




0.201 


-90.000 


-21.124 


o.c 


110. 




C.2C2 


-90.000 


-21.084 


0.0 


1 15. 




0.215 


-90.000 


-20.528 


0.0 


120. 




0.238 


-90.000 


-19.674 


o.c 


125. 




0.265 


-90.000 


-18.715 


o.c 


130. 




0.296 


-90.000 


-17.778 


0.0 


135. 




0.326 


-90.000 


-16.933 


o.c 


lAC . 




0.35A 


-90.000 


-16.208 


c.o 


1A5. 




0.379 


-90.000 


-15.610 


o.c 


15C. 




C.AOl 


-90.000 


-15. 132 


0.0 


155. 




C.A18 


-90.000 


-14.760 


c.o 


160. 




0. A32 


-90.000 


-14.482 


o.c 


165. 




C.AA2 


-90.000 


-14.282 


0.0 


170. 




C.4A9 


-90.000 


-14.148 


o.c 


175. 




0.A53 


-90.000 


-14.072 


G . 0 


180. 




C.A5A 


-90.000 


-14.048 


DIRECTIVITY/ISO 


F I 


3DB 


ELEV 3DB 





S.820CB A3. 25. 

THIS CIRECnVITY FIGURE V^AS COMPUTED 

ASSUMING THAT THE PATTERN IS WELL-BEHAVED. 
IT HAS NO MEANING IF THERE ARE SIDE-LOBES. 
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ALL LMTS ARE RKS 
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100 



OVERALL LENGTH OF STRICTURE 0.5A1 



CURRENT CISfRIHUTICN 1 FREQ= 0.230E 10 



ELEV 


F I 


E 


PSI 


EDB 


0.0 


0.0 


2.A70 


90.000 


-0.000 


10. 


r,C 


2.A00 


90.000 


-0.250 


20. 


0.0 


2.201 


90.000 


-1.002 


30. 


0.0 


1.9C3 


9C.000 


-2.266 


^0. 


0.0 


1. 5A7 


90.000 


-A.06A 


bC. 


O.C 


1.178 


90.000 


-6. A35 


60. 


c.c 


0.831 


90.000 


-9. A67 


70. 


0.0 
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THIS DIRECTIVITY FIGURE V«iAS COMPUTED 
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THIS DIRECTIvrTY FIGURE WAS COMPUTEn 

ASSUMING THAT THE PATTERN IS WELL-BEHAVED. 
IT HAS NC WEANING IF THERE ARE SIDE-LOBES. 
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II, Program SMIGUEL 



This program also is written in code FORTRAN IV for the IBM 360 
system. 

The program uses seven subroutines, as explained in SAMPLE PROGRAM 
SMIGUEL that follows. The subroutines supply to the main program the 
characteristics of the line and the non-uniformity corrections. The main 
program uses equations 33, first to find the input impedance, and then to 
find the current distribution. 

For each set of input data the impedance and current distribution are 
computed twice — one time without accounting for mutual impedance effects, 
and a second time accounting for them. 

Again, the program has a large number of comment cards that, hope- 
fully, make it self-explanatory. 
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phase, and very few models have been fully analyzed. The present paper is an 
attempt to analyze the operation of a zig-zag model that has the property of 
being symmetrical, and suitable for operation against ground. The radiation 
pattern of the antenna is obtained for different models of current distribution, 
and, finally, the impedance characteristics and an approximate current distribu- 
tion are obtained, using non-uniform transmission line theory. The results 
obtained show reasonable agreement with experimental data, and confirm 
conclusions drawn from physical considerations. 
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